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/A/'ms. We attempt to detect starlight reflected from the hot Jupiter orbiting the main-sequence star r Boo, in order to determine the 
albedo of the planetary atmosphere, the orbital inclination of the planetary system and the exact mass of the planetary companion. 
Methods. We analyze high-precision, high-resolution spectra, collected over two half nights using UVES at the VLT/UT2, by way of 
data synthesis. We interpret our data using two different atmospheric models for hot Jupiters. 

Results. Although a weak candidate signal appears near the most probable radial velocity amplitude, its statistical significance is 
insufficient for us to claim a detection. However, this feature agrees very well with a completely independently obtained result by 
another research group, which searched for reflected light from r Boo b. As a consequence of the non-detection of reflected light, 
we place upper limits to the planet-to-star flux ratio at the 99.9% significance level. For the most probable orbital inclination around 
i = 46°, we can limit the relative reflected radiation to be less than e = 5.7 x 10~ 5 for grey albedo. This implies a geometric albedo 
smaller than 0.40, assuming a planetary radius of 1.2 Rj up . 

Key words. Methods: data analysis - Techniques: radial velocities - Stars: individual: r Boo - planetary systems 



1. Introduction 

Since the detection of the first exoplanet orbiting a solar-type 
star, more than 400 exoplanets have been detected. The exis- 
tence of most of these planets was established by monitoring ra- 
dial velocity (RV) variations of the host star, originating from the 
gravitational pull of the unseen planetary companion. So-called 
hot Jupiters are giant planets only a few solar radii away from 
their host stars that provide the opportunity to attempt the detec- 
tion of starlight reflected from these planets. Five extended cam- 
paigns for the search for refected light by using high-resolution 
spectroscopy were completed by different groups (Charbonneau 
et al. 1999, Collier Cameron et al. 1999, Collier Cameron et al. 
2002, Leigh et al. 2003a., Leigh et al. 2003b, Rodler et al. 2008). 
Apart from Collier Cameron et al. (1999), who claimed a detec- 
tion of reflected starlight, which was later withdrawn (Collier 
Cameron et al. 2000), all campaigns resulted in a non-detection 
of reflected starlight, and upper limits to the planet-to-star flux 
ratio and to the geometric albedo of these planets were estab- 
lished: to date, the tightest 99.9% confidence upper limits on the 
geoemtric albedos of the hot Jupiters r Boo b, HD 75289b and 
v And b are 0.39 (Leigh et al. 2003a), 0.46 (Rodler et al. 2008), 
and 0.42 (Collier Cameron et al. 2002), respectively. These re- 
sults provided important constraints on models of the planetary 
atmospheres such as those by Sudarsky et al. (2000, 2003). As a 
result, models that predicted a high reflectivity for the planetary 
atmosphere could be ruled out for some of the studied planets. 

More recently, the albedos of several transiting hot Jupiters 
at visual wavelengths could be further constrained from mea- 
surements of the secondary transit events. For the transiting hot 
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Jupiter HD 209458b, Rowe et al. (2008) placed a very stringent 
3cr upper limit on the geometric albedo of 0. 17 in the wavelength 
regime of 400-700 nm. Using data of the CoRot-satellite of the 
transiting hot Jupiter CoRot-lb, Snellen et al. (2009) measured 
a phase-dependency of the planetary flux and reported an up- 
per limit on the geometric albedo of 0.2 in the wavelength range 
of 400-1000 nm. This upper limit of the geometric albedo of 
CoRot-lb was confirmed by an independent analysis by Alonso 
et al. (2009a). Using the same data set for the same target, 
Rogers et al. (2009) reports an estimate of the planetary geo- 
metric albedo to be 0.05. Alonso et al. (2009b) placed a very 
stringent upper limit on the geometric albedo of 0.06 of the hot 
Jupiter CoRot-2b in the wavelength range 400-1000 nm. The 
analysis of measurements of the secondary eclipse of the hot 
Jupiter HAT-P-7b with the EPOXI spacecraft and the Kepler- 
satellite led to an estimate of the geoemtric albedo of 0.13 at 650 
nm (Christiansen et al. 2009). Measurements of the secondary 
transit of the planet Ogle-Tr-56 in the z'-band indicate a low ge- 
ometric albedo less than 0.15 (Sing & Morales 2009). For a gen- 
eral summary of the secondary transit measurements we refer to 
Cowan & Agol (2010). 

t Bootis (HD 120 136 A) is a main-sequence star of spectral 
type F7, located at a distance of 15.6 pc from our Solar System. 
Butler et al. (1997) detected a massive hot Jupiter orbiting r Boo 
via RV measurements; we note that this star is one of the bright- 
est stars in the sky harboring a planet. Shortly after the discovery 
of the planetary companion, two research groups started cam- 
paigns for the search for reflected light by way of high-precision 
spectroscopy, which finally resulted in non-detections. The tight- 
est published 99.9% confidence upper limit to the geometrical 
albedo is p < 0.39 under the assumption that the planetary ra- 
dius R p = 1.2 Rj np , orbital inclination values i > 36° and grey 
albedo (Leigh et al. 2003a). These authors report the detection 
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of a candidate signal of marginal significance with a projected 
orbital RV semi-amplitude K p = 97 km s , which corresponds 
to an orbital inclination i « 40°. 

Here, we present our analysis of new observations of the 
planetary system of t Boo taken with the UV- Visual Echelle 
Spectrograph (UVES) mounted on the VLT/UT2 at Cerro 
Paranal in Chile. Section 2 describes the basic ideas of the search 
for reflected light. Section 3 provides an overview over the plan- 
etary system, while Section 4 outlines the acquisition and reduc- 
tion of the high-resolution spectra. Section 5 provides a detailed 
description of the data analysis with a data modeling approach. 
Finally, in Section 6 we present the results, which are discussed 
in Section 7. 



where 

««) = 0.09 + 2.39 - 0.65 (3) 

and the phase angle a is given by 

cos a = - sin i cos Imp. (4) 

Finally, the flux of the reflected light from the planet at the orbital 
phase <p follows from Equations Q] and [2] 

f( t />,i,A) = e(A)fi(cf>,i). (5) 



2. Starlight reflected from the planet 

High resolution spectroscopic searches utilize the fact that the 
spectrum reflected from the planet is essentially a copy of the 
rich stellar absorption line spectrum except for the following dif- 
ferences: 

- It is scaled down in intensity by more than five orders of 
magnitude in the visual for hot Jupiters (Section l2~TT i. 

- It is shifted in wavelength according to the relative orbital 
radial velocity of the planet (Section l2~2l . 

- It displays a different degree of rotational broadening corre- 
sponding to the rotational velocity of the star as seen from 
the hot Jupiter whose own rotation typically contributes very 
little to the line broadening (Section l2~3l 



2.1. Photometric Variations 

For exoplanets, the enormous brightness contrast between the 
star and the planet constitutes a considerable challenge when at- 
tempting to observe some kind of direct signal from the planet. 
For close-in planets, such as hot Jupiters, the main contribution 
to the flux in the visual consists of the reflected starlight and not 
the intrinsic luminosity (Seager & Sasselov 1998). 

According to Charbonneau et al. (1999), the amount of 
starlight reflected from a planet which is fully illuminated can 
be described by 



e(.l) = />(.!) | ^ 
' a 



(1) 



where p(A) denotes the albedo of the planet as a function of the 
wavelength A, R p the planetary radius and a the star-planet sepa- 
ration. The planetary radius of r Boo b is unknown; it can, how- 
ever, be estimated from a comparison with the transiting plan- 
ets, which provide an exact determination of their masses and 
radii (Bakos et al. 2007). The orbital radius is tightly constrained 
through Kepler's third law. 

As the planet orbits its star, the fraction of the illumi- 
nated part of the planet changes relatively to the observer. 
Consequently, the observed reflected light is reduced, depending 
on the model describing the scattering behaviour of the atmo- 
sphere, its orbital inclination i e [0°,90°] and the orbital phase 
<p e [0, 1] of the planet. We note that we adopt the convention 
that = represents inferior conjunction of the planet (for 
i = 90°, it would be the transit position). Lacking observational 
data for hot Jupiters we apply an empirical scattering model of 
the atmospheres of Jupiter and Venus (Hilton 1992), which can 
be approximated by 



fi(<f>,i) = l0- 0Ai{a \ 



2.2. Doppler shifts 

The planet orbiting its host star does not only produce a flux vari- 
ation (Equation |5j, but also a Doppler shift of the stellar spec- 
trum reflected from the planet. The RV semi-amplitude K p of 
that shift depends on the orbital inclination i, which is unknown 
for non-transiting planets. K p can be expressed by 



M s 



M p sin/ 



sin i, 



(6) 



where K s is the RV semi-amplitude of the star, and M s and 
M p sin i are the stellar mass and the minimum mass of the planet, 
respectively. The largest possible amplitude /r p , ma x occurs at or- 
bital inclination i - 90° (Figure[TJ. 

The instantaneous RV shift of the planetary signal with re- 
spect to the star depends on the orbital phase <p, 



V p = K v sin 2n(j). 



2.3. Rotational broadening 



(7) 



As the star rotates, each absorption line is subjected to Doppler 
broadening since individual points on the visible disk of the star 
have different instantaneous radial velocity. Equation[8]describes 
the projected rotational velocity of the star: 



, * = v* sin i - 2jt sin i 



(8) 



where R* and P 10t .+ are the stellar radius and the rotation period 
of the star, respectively. If the inclination of the rotational axis is 
unknown, the true velocity v of the star remains unknown. 

Mathematically, the rotational broadening constitutes a con- 
volution of the intrinsic stellar line profile with a half ellipse 
whose width is equal to ±v sin i (Gray 1992). The v sin i of the 
reflected stellar absorption lines, as they could be observed from 
the planet, can be calculated by way of Equation[9] 



"★,p 



2ttR+ 



1 



Pmt,* 



1 



orb 



(9) 



(2) 



where P 01 -b is the orbital period of the planet. Equation|9]is based 
on the assumption that the orbital plane and the stellar equa- 
tor are co-aligned. Analyses of Rossiter-McLaughlin effect RV- 
curves (Rossiter 1924; McLaughlin 1924) in transiting planets 
have shown that such a spin-orbit alignment is the case for the 
majority of the monitored systems (Fabricky & Winn 2009), al- 
though such systems have been recently discovered (e.g. Pont et 
al. 2009, Narita et al. 2009). 
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Fig. 1. RV curves of the planet orbiting t Boo for selected or- 
bital inclinations i. The maximum possible RV-shift of K p = 
156 km s _1 occurs at i = 90°. For lower inclinations, the mass of 
the companion has to be larger to produce the stellar reflex mo- 
tion with an RV semi-amplitude K s = 0.461 1 km s _I . The most 
probable orbital inclination is i = 46° under the assumption of 
a tidally locked system. The orbital inclination of i = 82° is 
the upper limit, coming from the non-transit detection, whereas 
i = 30° is the lower limit. For i < 17.2°, the companion would 
exceed the mass of 13 Mj up and therefore be a brown dwarf. 



To calculate the broadening coming from the planetary rota- 
tion, we use EquationfTOl 



Voj, P = v p sm 1 = 2n sm 1 „ — > 

* iot,p 



(10) 



where R p is the planetary radius and P ro t,p the planetary rotation 
period. We note that this equation is based on the assumptions 
that the planetary rotation axis is co-aligned with the orbital in- 
clination axis, which is furthermore co-aligned with the stellar 
rotational axis. We justify these assumptions by the following 
argument: Close-in planets such as hot Jupiter seem to be tidally 
locked with their parent stars (c.f. Shkolnik et al. 2005; Knutson 
et al. 2007). A tidal lock hypothesis suggest that the orbital incli- 
nation is co-aligned with the planetary rotation axis and with the 
stellar rotation axis. The latter assumption was confirmed by the 
aforementioned spin-orbit measurements of transiting planets. 
The two contributions v* >p (Equ. [9| and v pro j iP (Equ. ITOb to 



the broadening are finally summed up to 



Vrefl 



„ + V ■ 
.P proj.p 



(ID 



3. r Boo and its planet 

Table Q] summarizes the parameters of the planet and its host 
star. We note in passing that the system also contains a faint M- 
dwarf component at a separatation of ss 230 AU from the pri- 
mary (Patience et al. 2002; Eggenberger et al. 2003). 



3.1. Radial velocity amplitude and orbital inclination 

Knowing the stellar mass M*, the planetary minimum mass 
M p sin i, and the RV semi-amplitude of the reflex motion of the 
star K+, the RV semi-amplitudes of the planet can be constrained 



Table 1. Parameters of the star t Boo and its planetary com- 
panion. Abbreviations for the references are: B06 = Butler 
et al. (2006) and references therein, B97 = Baliunas et al. 
(1997), G98 = Gonzalez (1998), H00 = Henry et al. (2000), 
L89 = Latham et al (1989), L03 = Leigh et al. (2003a) , 
VF05 = Valenty & Fischer (2005). 



Parameter 


Value 


Error 


Ref. 


Star: 








Spectral type 


F7 IV-V 




B97 


V (mag) 


6.35 




VF05 


rf(pc) 


15.60 


0.17 


VF05 


M+ (Me) 


1.33 


0.11 


VF05 


R* (Ro) 


1.31 


0.06 


G98 


Pro. (d) 


3.31 




B97 


v sin; (km s~') 


14.9 


0.5 


H00 


Age (Gyr) 


1.3 


0.4 


VF05 


Planet: 








M p sini(M Jup ) 


4.13 


0.34 


B06 


fl(AU) 


0.0481 


0.0028 


B06 


e 


0.023 


0.015 


B06 


K* (kms- 1 ) 


0.4611 


0.0076 


B06 


Orbital period (d) 


3.31245 


0.00003 


L03 


T^o (ID) 


2 451 653.968 


0.015 


L03 



to ranging from K p = to 155.6 + 18.3 kms" 1 (Equation |6j 
Figure [TJ. 

Due to the absence of transits in high-precision photometry 
(Henry et al. 2000), we can constrain the range of possible or- 
bital inclinations. The minimum orbital inclination i that a transit 
event occurs can be calculated by 



R* + R P 



(12) 



where and R p is the stellar and the planetary radius, respec- 
tively, and a the semi-major axis. Using the parameters listed in 
Table 1, and the assumption that R p = 1.2 Rj up , we can exclude 
orbital inclinations of i > 82° for r Boo b. 

Baliunas et al. (1997) found that the star r Boo rotates 
rapidly with a period commensurate with the orbital period of 
the planet, suggesting tidal locking. This hypothesis seems to 
be very likely for close-in planets (e.g. Shkolnik et al. 2005; 
Knutson et al. 2007), but has not been observationally confirmed 
so far. A tidal lock enables us to place an estimate on the orbital 
inclination (Equation [T3l) under the assumption that the stellar 
equator and the orbital plane are co-aligned. 

In order to calculate the orbital inclination, we need to solve 



sin i 



•'proj,* 



2nR* 



(13) 



where v* is the measured rotational velocity of r Boo (v pro j ; * = 
v*sin/ = 14.9 kms 1 ) for the unknown orbital inclination i, 
and v the true, but unknown rotational velocity. This rotational 
velocity v+ can be estimated by way of Equation [8] which re- 
quires the knowledge of the stellar radius and stellar rotation pe- 
riod. Adopting the parameters listed in Table [1] we determine 
the maximum rotational velocity of the star to be v* = 20.0 + 
10 km s _1 . This gives an orbital inclination of i = 46^g °. The 
resulting mass of the companion to t Boo of M p = 5.7*jj Mj up 
would clearly confirm the planet hypothesis. We note that the er- 
ror values of i and M p correspond to the maximum error range. 
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3.2. Rotational broadening 

We investigated how the reflected spectrum from the planet 
would be affected by a tidal lock and planetary rotation, using 
the following asumptions: The hot Jupiter orbiting t Boo rotates 
in 3.1 days (1:1 resonance with the orbital motion), and the 
rotation axis is co-aligned with that of the star. In addition, the 
rotation of the planet is assumed to be prograde, and the plane- 
tary radius is R p = 1.2 Rj up . 

A tidal lock causes that an imaginary observer sitting on the 
planets always sees the same side of the star; the star appears 
to be non-rotating. For this reason, v+, p = km s _1 . According 
to Equation [TUJ we find that the reflected absorption lines from 
the planet are broadened by approximately v p sin i = 2 km s _1 . 
In order to estimate the resulting full width at half maximum 
(FWHM) of the reflected stellar absorption lines, we investigate 
slowly rotating stars of similar spectral type. 

A very good candidate is the star HD 136351 (spectral 
type: F6 IV), which shows a rotational broadening of v* sin i = 
2.8 km s _I (Reiners & Schmidt 2003). The average FWHM of 
the absorption lines is 12 kms -1 . Since the rotational velocity 
of the planet is estimated to be smaller (see above), we estimate 
the resulting FWHM of the reflected stellar absorption lines to 
be about 11.8 km s _1 . 



4. Data acquisition and data reduction 

We observed r Boo during two consecutive nights in June 2007 
using UVES (Dekker et al. 2000) mounted on the VLT/UT2 
and collected a total of 406 high-reolution spectra (Table [2]). 
In addition to that, we took spectra of the slowly rotating 
star HD 136351 and the rapidly rotating B-star HD 116087 
(v* sin / = 241 km s _1 ; Hoffleit & Jaschek 1991). The latter star 
was observed for the identification of telluric features in the red 
part of the visual. The dates of the observations were selected 
in such a way that the observations were carried out at orbital 
phases at which the planetary signal was strongest, i.e. close to 
the position of superior conjunction, which would be the sec- 
ondary transit in case of i = 90°. We aimed at taking observa- 
tions in the phase ranges (p = 0.30 to 0.45 and 0.55 to 0.70, but 
avoiding (f> = 0.50 because of intense blending of the planetary 
and stellar absorption lines. 

The observations were conducted using the red arm of UVES 
with a non-standard setting centered at wavelength A = 430 nm 
providing 47 full orders and covering the wavelength range 
A = 425 to 632 nm. We observed through a 0.3" slit and the im- 
age slicer IS#3, providing us a spectral resolution of R = A/AA = 
110 000. The integration time for our target was selected to pro- 
vide the maximum count rate but at the same time avoiding to 
reach the non-linear regime of neither one of the two CCDs. The 
exposure times for our bright target ranged between 60 and 400 s 
(120 s on average) and 30 and 60 s (40 s on average), respectively 
for the first and the second night. We carried out our observations 
in the fast-R/O mode; the total dead-time between two exposures 
was 16 s. Due to clouds and thick cirrus, we only collected half 
of the expected spectra. 

In addition to the science frames, we obtained a large num- 
ber of calibration exposures. Before the start of the first night, we 
took 163 flat-field and 104 bias exposures. In the afternoon be- 
fore the second night, we collected 245 flat-field and 65 bias ex- 
posures. The large number of calibration frames was important 
to achieve a low photon noise in the finally combined flat-field 
frame (masterflat). 



4.1. Data reduction 

In order to prevent introducing data reduction artefacts to the 
data, we aimed at keeping the spectra as raw as possible. 
Consequently, we only adopted the absolutely necessary data re- 
duction steps, as pointed out in the following. Effects like er- 
rors in the wavelength calibration, trends in the continuum, in- 
strumental profile changes were considered in the model for the 
star/planet (see Section 5). 

In the first step, we created masterbias frames for each night, 
being the median of all bias exposures per night. These mas- 
terbias frames were then subtracted from the science frames as 
well as from the flat-field frames. The error of the flux in the 
science frames was determined on the basis of Poisson statistics 
and propagated in the further data analysis steps. In the follow- 
ing step, for each night the flat-field exposures were scaled with 
their inverse exposure times and then combined into their me- 
dian (masterflat) to permit high-quality flat-field correction in 
order to compensate for sensitivity variations of the pixels. The 
science frames were divided by the appropriate masterflat. Since 
both the flat-field frames as well as the object frames show a sim- 
ilar Blaze function, this function was then mostly removed from 
the object frame. Finally, 1 -dimensional spectral orders (pixel 
vs. flux) were extracted from the 2-dimensional frames. 

We retrieved 47 orders of 4096 pixels each. No order merg- 
ing was applied. The observation of the UVES Thorium-Argon 
(Th-Ar) lamp enabled us to assign each pixel the appropriate 
value of the wavelength. For each night, we used the Th-Ar spec- 
trum observed before/after the science exposures and established 
an 8th order polynomial dispersion solution. All these data re- 
duction steps were performed with the MIDAS software pack- 
age. 

5. Data analysis 

After the data reduction, we identified cosmic-ray hits by way 
of the following procedure: for each spectrum, we compared the 
flux in every pixel with the median flux of the same pixel in 
the three predecessor and the three successor spectra, which had 
been scaled to the same flux as the spectrum under considera- 
tion. We flagged those pixels where the difference exceeded 6<x 
as cosmic-ray hits. These pixels were then excluded from further 
analysis. We furthermore discarded the most weakly exposed re- 
gions of each echelle order (the first 400 pixels as well as the 
last 240 pixels) since the flux level of the continuum was sub- 
jected to strong variations. In addition, the spectral order around 
the wavelength of A — 630 nm was excluded from further anal- 
ysis, since data were affected by the strong telluric OI features. 
In order to speed up the data analysis, we co-added spectra into 
sufficiently narrow phase bins in such a way that phase smearing 
of the stellar absorption lines originating from stellar RV varia- 
tion, sub-pixel shifts and the effect of barycentric motion of the 
Earth remained negligible. The main criterion for the size of the 
phase bins was that the unseen planetary lines did not suffer from 
smearing in excess of 2 km s _1 . With that step, we were able to 
reduce the number of spectra from 406 to 58. 

5.1. Data synthesis method 

We model the starlight reflected from the planet as a copy of 
the stellar spectrum, strongly scaled down in brightness and 
Doppler-shifted according to the orbital motion of the planet. 
The target spectra have an average S/N of 300 to 600 per spec- 
tral bin. With expected planet-to-star flux ratios of the order of 
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Table 2. Journal of observations. The UTC times and the orbital phases of the planet are shown. Ephemerides used are: 
BJD 2451653.968 + 3.31245 x E (Leigh et al. 2003a quoting Marcy priv. comm.). 



Night 


UTC start 


BJD start 


<t> 


UTC end 


BJD end 





spectra 


1 

2 


2007/06/16 23:24 
2007/06/17 22:45 


2454268.4773 
2454269.4498 


0.30 
0.59 


2007/06/17 04:12 
2007/06/18 04:26 


2454268.6784 
2454269.6873 


0.36 
0.66 


100 
306 



a few times 10~ 5 , it is clear that the reflected spectrum from the 
planet is deeply buried in the noise of the stellar spectrum. The 
weak planetary signal is boosted by the large number of spectra, 
and more importantly, by the combination of the approximately 
1500 absorption lines, achieved using the data-synthesis method 
(cf. Charbonneau et al. 1999) described below in a nutshell (for 
a detailed description see Rodler et al. 2008). 

In the first step, a high S/N, virtually planet-free superspec- 
trum is computed by co-adding the observed spectra after cor- 
recting their wavelength values for the barycentric motion of 
the Earth and for the stellar orbital motion. Then, we construct 
a model to describe each of the original, unmodified object 
spectra: the dominant stellar signal is represented by the su- 
perspectrum, shifted to the position of the observed spectrum. 
Imperfections in flux and wavelength are corrected by using the 
approach described in Rodler et al. (2008). Furthermore, the con- 
tribution to the spectrum of the planetary signal is created as 
follows: we adopt the spectrum of the slowly rotating F6 IV 
star HD 136351 to mimic the expected sharp reflected lines of 
the planet (Section [3~2l . For each observed object spectrum, we 
need to co-align that spectrum of HD 136351 with the stellar 
model spectrum (shifted superspectrum). The co-aligned spec- 
trum of HD 1 3635 1 is then scaled down by the factor e(A) fi(<f>, i) 
and shifted by velocity V p (K p , <p) with respect to the stellar spec- 
trum. Hence, the two free parameters of the planetary model are 
the planet-to-star flux ratio for the fully-illuminated planet e(A), 
and the orbital inclination ;, which corresponds to an RV semi- 
amplitude of the planet K p - A" p>max sin! = 155.6 sin; kms -1 . 
Concerning the albedo function p(A), we adopt (i) a grey albedo 
assumption (i.e. p(A) = constant for all wavelengths A) and 
(ii) the irradiated atmospheric Class IV model by Sudarsky et 
al. (2000), which predicts higher reflectivity at shorter wave- 
lengths (Fig.|2]i. 

The search range for the RV semi-amplitude comprises K p = 
40 to 180 km s _1 (corresponding to orbital inclinations i = 15° 
to 90°, plus twice the error of K p max (i.e. 18.3 km s _1 with a 
step width of 3 kms -1 . This is a good compromise between 
computing time and sampling the expected average reflected ab- 
sorption line profile with the FWHM of a: 12 km s -1 . Since the 
tidal lock hypothesis and a spin-orbit alignment is very likely, 
low-inclination orbits can be furthermore excluded: the observed 
vsin/ = 15 kms -1 would imply an improbable true rotational 
velocity above 50 km s for an F7 IV-V star at very low or- 
bital inclinations (cf. rotational velocities of late F-type stars; 
Glebocki & Stawikowski 2000). Using simulations, we found 
that for small inclinations of the planetary orbit, where the plan- 
ets appear only slightly illuminated, and the method would fail 
to detect Jupiter-size objects even with very high albedos. 

For the second parameter e(Ao), we search for the planetary 
signal in the range 10~ 4 to 10~ 5 with a step size of 5 x 10~ 6 . 

5.2. Determination of the confidence level 

Once the best model M[K p , £{A)] has been evaluated, we deter- 
mine the confidence level of the x 2 minimum. We do not infer 
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400 450 500 550 600 

X [nm] 

Fig. 2. Different albedo spectra of atmospheric models (taken 
from Sudarsky et al. 2000) are shown. The irradiated (dots) and 
isolated (dashed) Class IV models describe atmospheres of plan- 
ets with temperatures T e $ « 1300 K. Contrary to the isolated 
model, the irradiated model assumes that no reflective clouds ex- 
ist in the upper layer of the planetary atmosphere, which results 
in a very low albedo. The solid line depicts a Class V roaster, 
describing the atmosphere of a planet with r e ff > 1500 K, hav- 
ing a high reflective silicate cloud deck in the upper layers of the 
atmosphere. This atmospheric model was already ruled out for 
t Boo b by Leigh et al. (2003a). 

the confidence from the probability of the minimum x 2 value, 
because due to unknown systematic errors coming from the tele- 
scope and instrument the true errors are largely underestimated 
when considering only photon statistics and detector read-out 
noise. Therefore no reasonable probability estimate is possible 
from the x 2 value alone. For this reason, we apply the bootstrap 
randomization method (e.g. Kiirster et al. 1997). Retaining the 
orbital phases, we randomly redistribute the observed spectra 
amongst the phases, thereby creating a large number (N = 3000) 
of different data sets. Any signal present in the original data is 
now scrambled in these artificial data sets. For all these ran- 
domised data sets, we again evaluate the model for the two free 
parameters, and locate the best fit with its specific x 2 minimum. 
We set m to be the number of best-fit models to the N randomised 
data sets that have a minimum^ 2 less or equal than the minimum 
X 2 found for the original data set. The confidence level can then 
be estimated by 1 - m/N. 

6. Results 

6.1. Grey albedo model 

Our data analysis adopting a grey albedo model resulted in a 
^-minimum at a planet-to-star flux ratio e = 3.3 x 10~ 5 and 
an RV semi-amplitude K p = 103 km s corresponding to an 
orbital inclination of i — 41°. Figure [3] (upper panel) shows the 
X 2 contour map in the e vs. K p plane. However, using bootstrap 
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Fig. 3. Contour maps of x 2 for the model parameters of the 
planet of K p (RV semi-amplitude of the planet) and e(/tn) (planet- 
to-star flux ratio); this plot shows the result of the analysis of the 
r Boo b data for two different atmospheric models. The contour 
levels follow the sequence C* = ^mi n + 10~ 6 0^ 2 -;^ in ) 3 for A: > ^> 
where y 2 . is the minimum y 2 . 

A mm A - 

Upper panel: grey albedo model (model #1). Lower panel: irra- 
diated Class IV albedo function (model #2). For both cases, a 
Venus-like phase function was adopted. The minimum x 2 was 
found at K p = 103 km s" 1 and e = 3.3 x 1(T 5 , and K p = 
103 km s _1 and e(/lo) = 2.5 x 10~ 5 , respectively, for model #1 
and model #2. A bootstrap analysis with 3000 trials revealed 
that these features are uncertain with a false alarm probability 
of 3.6% and 7.3%, respectively, for model #1 and model #2. For 
comparison, we show the candidate feature found by Leigh et 
al. (2003a) at ^p = 97 km s _1 and e = 4 x 10~ 5 , assuming the 
grey albedo case (upper panel, cross). This candidate feature has 
a FAP of 1.4%. Note that the corresponding geometric albedo p 
is shown under the assumption that R p — 1.2 Rj up (the right-hand 
y-axis). 



randomisation with 3000 trial data sets we find this ^-minimum 
to be uncertain with a false alarm probability (FAP) of 3.6%, and 
consequently do not consider it as a detection of reflected light 
from the planet. 
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Fig. 4. Contour map showing confidence levels for the upper lim- 
its to the planet-to-star flux ratio e(/tn) as a function of to the RV 
semi-amplitude K p of the planet (lower x-axis), or orbital in- 
clination i (upper x-axis). From top to bottom, these levels are: 
99.9, 99.0 and 90.0 % confidence. For the planetary model we 
assumed a Venus-like phase function as well as a irradiated Class 
IV albedo function (dashed lines) and a grey albedo (solid lines). 
The corresponding geometric albedo p is shown under the as- 
sumption that R p — 1 .2 Rj up (the right-hand y-axis). 



6.2. Irradiated Class IV model function 

The analysis with the irradiated Class IV albedo model did not 
yield any evidence for reflected light either. We find the formal 
^-minimum at a planet-to-star flux ratio e(/lo) = 2.5 x 10~ 5 
and an RV semi-amplitude of K p = 103 kms , which again 
corresponds to an orbital inclination i = 41°. Figure [3] (lower 
panel) shows the x 2 contour map in the e vs. K p plane. From 
bootstrap randomisation with 3000 trials we found this feature 
to be uncertain with a FAP of 7.3%. 



6.3. Upper limits 

Given the non-detection of the planetary signal, we determined 
upper limits of the planet-to-star flux ratio e(Ao) for different 
confidence levels as a function of the RV semi-amplitude K p . 
Figure|4]shows that the upper limits to the planet-to-star flux ra- 
tio decrease with increasing orbital inclination, which is a direct 
consequence of the illumination geometry. As can be seen from 
Figure |4] (solid lines), we have the highest sensitivity for detec- 
tion of the planetary signal at high orbital inclinations. For the 
most probable orbital inclination of 46°, the 99.9 % confidence 
upper limit to the planet-to-star flux ratio for the grey albedo 
model is e = 5.7 x 10~ 5 , while it is e = 6.5 x 10~ 5 for the ir- 
radiated Class IV albedo function. For an orbital inclination of 
i = 60°, the 99.9 % confidence upper limit to the planet-to-star 
flux ratio for the grey albedo model is already e = 5.1 X 10" 5 , 
while it is e = 5.7 x 10~ 5 when adopting the irradiated Class 
IV albedo function. Similar to our results for the hot Jupiter 
HD 75289 Ab (Rodler et al. 2008), Figuregjshows that the upper 
limits established by adopting the grey albedo model are deeper 
than the ones found with the irradiated Class IV model. 

Assuming a planetary radius R p = 1.2 Rj up and an orbital 
inclination i a 46° (cf. Section |3), we find the upper limit to the 
geometric albedo to be p < 0.40 for the grey albedo model and 
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p < 0.44 for the irradiated Class IV model at a wavelength of 
Ao = 462 nm, which corresponds to the centre of gravity of the 
irradiated Class IV albedo function in the observed wavelength 
range. For comparison, at this wavelength Jupiter's geometric 
albedo is p = 0.44 (Karkoschka 1994). 

7. Discussion and conclusions 

We have observed the hot Jupiter orbiting t Boo for two half 
nights with UVES, mounted at the VLT/UT2, in an attempt to 
measure starlight reflected from the planetary companion. 

- The data analysis using the data synthesis method did not 
reveal a detection of the reflected light from the planet. 
The signal-to-noise ratio of the combined data was not 
sufficient to significantly detect the planetary signal; instead 
we placed upper limits to the planet-to-star flux ratio for 
different possible orbital inclinations and confidence levels. 
Under the assumption of a planetary radius of R p — 1 .2 Rj np , 
we rule out atmospheric models with geometrical albedos 
p > 0.40 for r Boo b. We confirm the finding by Leigh et al. 
(2003a) that the upper limit to the geometric albedo of the 
planet (e.g. p < 0.39 for orbital inclinations ; > 36° and grey 
albedo) implies that this planet has a lower reflectivity than 
the atmospheres of Jupiter (p = 0.44) and Saturn (p = 0.62). 

- For both albedo models, we find a formal ^-minimum at 
K p = 103 kms -1 . For the analysis adopting the irradi- 
ated Class IV albedo model, we determine the FAP of the 
candidate feature to be 7.3%. In comparison, Leigh et al. 
(2003a) finds a candidate feature of marginal significance 
(FAP=0.9%) with K p = 90 km s" 1 . 

For the analysis adopting the grey albedo model, we find that 
this candidate feature has a low FAP of 3.6% with a planet- 
to-star flux ratio of e = 3.3 x 10~ 5 . If genuine, this would 
correspond to an orbital inclination of i = 41° and a geo- 
metric albedo of p — 0.23 assuming R p = 1.2 Rj up . For the 
grey albedo model, Leigh et al. (2003a) reports about the de- 
tection of a non-significant candidate feature with a similar 
value of K p = 97 km s _1 , having a low FAP of 1.4%. 
These similar, but non-significant candidate features were 
found by using completely independent data sets and data 
analysis methods. Leigh et al. (2003a) analyzed high- 
resolution spectroscopic data, which had been obtained dur- 
ing 17 nights with the Utrecht Echelle Spectrograph (UES), 
mounted on the William Herschel-telescope in La Palma, 
Spain. Contrary to our analysis strategy, they used a decon- 
volution approach (Collier Cameron et al. 2002) for the ex- 
traction of the planetary signal. In this work, we analyzed 
UVES data, obtained during two half nights, by using the 
data synthesis method. These results are also consistent with 
an estimate of the orbital incliation of i as 46° (correspond- 
ing to a semi-amplitude of K p = 112 km s ) and assuming 
tidal locking. 

Leigh et al. (2003a) achieves a higher sensitivity by adopting 
the irradiated Class IV model atmosphere, whereas our high- 
est sensitivity is attained with the grey albedo assumption. 
One explanation for this might be that Leigh et al. (2003a) 
were using data with a slightly better wavelength coverage at 
blue wavelengths (407-647 nm), while our data set covered 
the wavelength range from 425 to 632 nm. However, it seems 
to be more likely that this discrepancy between the best-fit at- 
mospheric models is just coincidence, since both campaigns 
resulted only in non-significant candidate features. 
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